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intensity monitor. Several calibration trials and three 
over-water experiments were conducted. A value on the 
order of 5x10"^ meters was found for the extinction 
coefficient in clear weather. 
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ABSTRACT 



A method of measuring laser beam extinction in the 
atmosphere over the ocean has been devised and satis- 
factorily utilized. An attempt has been made to correlate 
the observed values with existing meteorological conditions. 

The apparatus consisted of equipment, already being utilized 

✓ 

for scintillation measurements, modified by the addition 
of a source intensity monitor. Several calibration trials 
and three over-water experiments were conducted. A value 
on the order of 5x10”^ meters was found for the 
extinction coefficient in clear weather. 
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I. PROBLEM DESCRIPTION 



A. GENERAL 

The increasing utilization of laser technology in long- 
range applications such as communications links , weapons 
designators /range finders, and high energy laser weapons 
makes it desirable to have the capability of predicting 
the effectiveness of the particular device at the termina- 
tion of the path through the atmosphere. A highly collimated 
beam of monochromatic light, such as laser light, is affected 
by many phenomena among which scintillation, beam wander, 
beam spread, and attenuation are only a few. The develop- 
ment of a procedure by which the latter, commonly called 
extinction, could be measured* in a path through the atmos- 
phere directly above the ocean surface was the goal of this 
experiment. 

Atmospheric attenuation (extinction) is generally 
caused by two phenomena: absorption by the gaseous mole- 

cules, and scattering by haze, fog, rain, and other aerosols. 
The absorption problem is minimized during system design 
by selecting a lasing medium that produces light of the 
proper wavelength so as to be within one of the several 
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transmission "windows*' in the atmosphere. Thus, the degree 
of extinction becomes a function of the meteorological 
conditions . 

Ideally, when sufficient data has been obtained, a 
system operator would only have to observe the existing 
meteorological conditions and thus predict the effective- 
ness of his equipment. Extinction measurements over 
terrain have previously been made /lief. l7. However, such 
measurements over the marine boundary are virtually non- 
existent . 

B. MATHEMATICAL DESCRIPTION 

If a beam of monochromatic flux enters an attenuating 
layer, the Bouger-Lambert Law /Ref. X] states that the 
fliix at any point within the layer is given by 

P= PoS'^’^ (1) 

where P= Flux remaining after the beam 
has traversed a path length % , 

Po = Flux entering the attenuating 
layer, 

= Extinction Coefficient, 

•pc = Path length in the attenuating 
layer. 

The flux described in equation (1) above has the units 
of energy per unit time over a given area and will be 
further described. 
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Considering both absorption and scattering the extinc- 
tion coefficient 

becomes = o< ^ (2) 

where <=*< = absorption coefficient, 

and ^ = scattering coefficient. 

In this experiment, 0.6328 micrometer wavelength laser 
light was utilized and since this is within one of the 
transmission "windows" no attempt was made to distinguish 
the relative values of *=< and ^ , but only to determine 
their combined value. 

The relationship (Equation 1) holds for a spherical 
wave. The detailed mathematical development, following 
references 3 and 4, presented in Appendix A, shows that a 
laser source producing pure monochromatic radiation can be 
considered to be a point source of spherical waves, pro- 
vided the observations are made in the paraxial region, at 
a distance from the source much greater than the off-axis 
position. 

Thus, considering again equation (1), and assuming that 
the effective source of the radiation is at the Gaussian 
waist of the diverged beam and radiates P® watts into a 
solid angle , then the intensity at a position near the 

beam axis at a distance R is given by 
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2 

where I = Intensity at R in Watts /meter 

2 

_Ti_ R = Area of beam at distance R. 
Consider propagation in a controlled atmosphere where 
scattering aerosols are at a minimum, so that , 

and a calibration factor can be determined from 

T = (4) 

r2 

where the subscript indicates all measurements taken during 
the controlled trial. Then from (4) 
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so 
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Taking the Natural Logarithm results in 

In (IR^) = In (I^Ro ^ - yU,R. (7) 

2 

A plot of In (IR ) versus R will yield the ordinate 
intercept 

In (I^R^ = In C. (8) 
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So in general for longer ranges where the determination 
of jul, is desired 
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and 



from which: 



In IR^ = In C y<.R 



In C - In IR^. 



( 10 ) 

( 11 ) 



Values of In /iR^J at each value of R can thus yield 
a calculated value of . 

Alternatively, it is only necessary to plot In £ IR^^ 
versus R for a series of measurements at increasing ranges 
and from this determine yU. as the absolute value of the 
slope. This assumes, of course, that yu. remains constant 
over the time and distance of the measurements. 

To obtain the necessary data for the plot of In / IR^J 
versus R the procedure briefly described below is followed: 

A horizontally fanned beam of constant intensity sweeps 
vertically across a detector located at a known distance R, 
4000 times each second. The peak value of each detected 
pulse is recorded and the pulse height distribution 
(spread out because of scintillation) is obtained. The 
average value of the pulse height is calculated and taken 
to be a measure of the average intensity I, at the distance 
R. 
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II. EXPERIMENTAL APPARATUS 



A. GENERAL 

Since the ultimate goal of the experiment was to measure 
the extinction coefficient over the marine boundary in con- 
junction with simultaneous meteorological measurements, the 
use of the Research Vessel ACANIA was necessary. Because 
of the limited availability of the ship's services, it was 
deemed desirable to test equipment and procedures thoroughly 
prior to actual shipboard measurements. This was accomplished 
first in the 120 meter long basement hallway of the Naval 
Postgraduate School's, Spanagel Hall. To accomplish this 
the laser source and scanning equipment (to be described 
in detail in the next section) were placed at one end of 
the hallway, and the detector and processing equipment were 
moved the length of the hallway in twenty meter increments. 

When operation of the equipment in this controlled 
environment was deemed satisfactory the apparatus was re- 
configured. The laser source and scanning equipment were 
mounted on a training and elevating mechanism in a two- 
wheeled trailer. The detecting and processing equipment 
was mounted in an instrumented bus (Fig. 1) which provided 
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Naval Postgraduate School Mobile Optical 
Research Laboratory. 




Interior, Showing 18'lnch Telescope. 



FIGURE 1. INSTRUMENTED BUS and 18-Inch Telescope 
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portability, shelter, and 110 V.A.C. power from a 3.5 K.W. 
gasoline engine driven generator. This configuration pro- 
vided the necessary support for extending the path length 
to 300 meters in 100 meter increments along O'HARE Avenue 
at the Naval Postgraduate School Annex. Operation in this 
manner was the final preliminary check prior to shipboard 
operation. For shipboard measurements the apparatus was 
again re-configured. The laser source and scanning equip- 
ment were mounted in the bus and the detecting and processing 
equipment were aboard ship. The reasoning behind this final 
re-configuration will become apparent in the following 
paragraphs . 

B. DETAILS OF EXTINCTION EQUIPMENT 

The laser source, scanning equipment, detection, and 
processing equipment utilized for extinction measurements 
were also utilized for scintillation and aperture averaging 
measurements; however, only that equipment in direct support 
of the extinction experiment will be described in detail. 

1. Bus Mounted Equipment 

Equipment described in the following paragraphs 
was mounted aboard the instrumented bus utilized by the 
optical propagation group. The bus itself was a converted 
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General Motors Corporation passenger bus. The forward 
one-third of the vehicle contained storage cabinets and 
workbench space and the drivers compartment. The after 
two -thirds of the vehicle contained electronic equipment 
racks, workbench, and an 18-inch telescope mounted on a 
platform capable of being raised clear of the floor on 
three hydraulic camper jacks extending through the floor 
to the ground. This arrangement provided a stable mount 
for the telescope and all equipment mounted on it and was 
thus isolated from any vibration caused by personnel moving 
about in the vehicle. A portable 3.5 kilowatt generator 
provided 110 V.A.C. power to all equipment, 
a. Laser Source 

The laser utilized as a source for the extinc- 
tion measurements was a 0.6328 micrometer He-Ne laser 
manufactured by C.W. Radiation INC. The output of the laser 
was fitted with a diverging lens controlled by a micrometer 
screw adjustment. The laser, scanning equipment, and source 
intensity monitor were mounted on a common plate attached 
to the right -hand- side of the bus mounted 18-inch telescope 
(Fig. 1) . This provided a means by which the laser could 
be trained, elevated, and aimed at the shipboard mounted 
detectors . 
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b. Scanning Equipment 

(1) Physical Description . Mounted directly 
in front of the laser was a mirror positioned at a 45° 
downward angle (Fig. 2) . This mirror performed two func- 
tions . First, it directed the beam down onto the scanning 
mirror. Secondly, it was capable of being bent, by means 
of a micrometer screw-adjustable knife-edge placed in the 
vertical axis, so as to produce a "FANNED” beam elongated 
in the horizontal plane. This arrangement was capable of 
producing a beam at the target ship sufficiently wide so as 
to provide uniform coverage of the detector as the ship 
moved and swung about its anchor. 

The actual scanner consisted of a small, 
1.0 cm X 1.0 cm, mirror mounted at a 45° upward angle 
^ef. 6, p. 4Q7 . This mirror, attached to two piano wires 
and mounted between the poles of a strong magnet, was set 
into oscillation and driven by a power amplifier. The 
amplifier controls were adjustable so as to provide a 
maximum vertical scan angle of 14 milliradians at a fre- 
quency of approximately 2000 HZ. The result of this beam 
processing was to allow detection of maximum intensity 
independent of the ship's vertical motion. 
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Cylindrical Mirror 
(micrometer adjustment provides 
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Retractable Mirror 

^ _ , (shown in position for 

FIGURE 2 Diagram of Beam Path sampling beam intensity) 
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(2) Mathematical Description . The following 



is a brief analysis of the scanned beam in terms of the 
vertical angular width of the beam and the scanning rate 
of the oscillating mirror. This permitted a prediction of 
the time width of the detected pulse which was then found 
to be in agreement with the measured value. To determine 
this the beam was fanned and both the unscanned and scanned 
beam sizes were measured at a convenient distance from the 
scanning mirror. Refer to Figure 3 for the details of the 
geometry; where ^ (the distance from the scan mirror) = 
2316.5 cm, h^ (the unscanned beam height at / ) =5.5 cm, 
and h .2 (the scanned beam height at / ) = 26 cm. 



of coverage, equals 2.37 milliradians and 2-0^, the angle of 
coverage of the scanned beam, equals 11.22 milliradians. 
Thus, -0^ =5.61 milliradians. 



From this c< , the unscanned beam angle 



Now the pulse width is given by 



where 



/V - ^ ^ 



( 12 ) 




scanning mirror by the 



detector. 



and ~ maximxim angular scan rate. 
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FIGURE 3 Geometry of Pulse Width Determination 
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for 0 where 






D equals the height of the detector. 



Thus , 



r- 



c< + 









(13) 



For large values of ^ , the second term 
in the numerator of (13) becomes negligible thus rendering 
the time width of the pulse nearly constant with range. 



Now 



-e- ^ 



5;o CD t) 



or 



d-e 



cj> -Q- Cos cD)t: . 



(14) 



(15) 



For maximum value of -O- , Cos uj't ~ dL 
and 



with 






(16) 



^ 3.000 ) 

and ■% =5.61 milliradians , as before; 

Thus, = 70.5 radians /sec. 

For a detector of two centimeter diameter 
and Jl as before, equation (13) yields a time width of the 
pulse, ^ , of 45.8 microseconds . This agreed with 

oscilloscope measurements of the detected pulses. 
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Figure 4 shows the scanning equipment and 
mounting details; the entire assembly was referred to as 
the "FAN-SCAN," 

c. Source Intensity Monitor 

Periodically throughout the extinction measure- 
ments it was necessary to sample the laser output as the 
beam left the scanning mirror to ensure that no intensity 
degradation was present at the source. To facilitate this, 
the Source Intensity Monitor was devised (Fig. 4). An 
adjustable mirror was installed which could be swung into 
position to interrupt the outgoing beam and direct it back 
to the monitor detector mounted above and behind the source 
laser. The detector, a PIN-5D Photodiode was mounted in 
an adjustable mount to permit fine alignment of the sensi- 
tive surface with the sampled beam. The output of the 
detector was displayed on a Monsanto OS 226 P, oscilloscope. 

2 . Shipboard Mounted Equipment 

The detection and signal processing equipment aboard 
the R.V. ACANIA included, detector, gyro-stabilized detector 
mount, PAR 113 Low-Noise Differential Preamplifier, Demodu- 
lator, Hewlett-Packard 7562A Log Voltmeter, Victoreen PIP- 
400 Multichannel Analyzer, and Hewlett-Packard HP 9810 
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Laser Source and Scanning Equipment Mounted 
on ISHnch Telescope. 




CloseMip Showing, from left, Monitor Detector, 
Laser, Cyl indrical Mirror, Oscillating Mirror, 
and Retractable Mirror. 



FIGURE 4."FAN-*SCAN and Source Intensity Monitor Arrangement 
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computer and associated plotter. The following paragraphs 
describe in detail the arrangement and fxinction of each of 
the above. 

a. Detector 

The detector utilized for the extinction measure- 
ments was a TIX L69 Silicon Avalanche Photodiode reversed 
biased to 159 volts. This provided an acceptable Signal- 
To-Noise Ratio at the extended ranges. The detector was 
mounted in a holder fitted with lens and an aperture adjust- 
able to 4.2 cm diameter. For more detailed discussion of 
the detector see ZRef. 77. 

b. Gyro Stabilizer 

The detector was mounted on the gyro stabilized 
mount. The stabilization was necessary to maintain align- 
ment of the detector with the laser source ashore. 

The gyro was an inertially stabilized platform 
utilizing two iron flywheels mounted on each end of a 1/2 hp 
capacitor start electric motor (Fig, 5). The shore based 
laser source utilized for the extinction measurements also 
served as the beacon onto which the gyro feedback circuit 
locked. Feedback was provided by a four-quadrant position 
and intensity photodetector with associated electronics. 
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Gyro unit;Showing top mounted Avalanche 
photodiode detector. 




Signal processing equipment; PARIIS’sCleftl; 
equipment in the rack ( top to bottom) ; Demodula- 
tor, Log Voltmeter/ Amplifier, Oscilloscope, and 
Multi-channel Analyzer. 



Hgure 5. Signal 



Detection and Processing Equipment 
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A small aiming telescope and "JOY STICK" control 
were mounted on the left side of the gyro to aid the gyro 
operator in aiming the control detector at the laser source. 
The gyro platform was positioned on the upper deck of ACANIA 
and was sheltered from the elements by an aluminum tool 
shed open on one side. The signal processing equipment 
described below was bench mounted in the "DRY LAB" on the 
main deck. (Fig. 5) . 

c. Preamplifier 

The output of the detector was fed through an 
F.E.T. follower (which served the function of impedance 
matching) to the PAR 113 Low-Noise Preamplifier. This 
provided high gain, 10 to 10^, low noise amplification, of 
the received signal. Adjustable High and Low Frequency 
Rolloffs allowed selection of the bandwidth. 

d. Demodulator 

The modulated signal was next passed to the 
Demodulator. This device sensed the scanning beam pulse, 
followed the pulse to its maximum and held that instanta- 
neous signal level in a flat-topped pulse until the next 
maximimi scanned across the detector at which time it held 
that maximum level, and so on. The output of this circuit 
was thus a series of step functions. 
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e. Log Voltmeter 

The Hewlett-Packard 7562-A Log Voltmeter was 
the next instrxmient in the processing step. At this point 
the Logarithm of the signal was obtained. 

f. Multichannel Analyzer 

The Victoreen, PIP-400, Multichannel Analyzer 
collected, stored for further readout, and displayed the 
sampled Log Voltage Pulses. Each of the available 200 
channels corresponded to a different pulse height and the 
ntjmber of pulses of any particular height received during 
the sampled period were stored in the respective channels. 

g. Computer and Plotter 

The output of the Multichannel Analyzer was 
fed into a Hewlett-Packard 9810A calculator. The calculator 
provided controlling functions to the associated plotter. 

The calculator sampled each channel of the Multichannel 
Analyzer, stored the information and plotted the pulse 
height distribution. The calculator then computed a least- 
squares polynomial fit of a Gaussian distributed curve to 
the data points. Using the position of this curve with 
respect to read-in calibration points and the Standard 
Deviation of the curve, the average intensity of the re- 
ceived signal was computed. See Appendix B for a detailed 
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analysis of the procedure for obtaining the average pulse 
height from the log voltage distribution. For a detailed 
description of the computer program see Ref. 8. Figure 6 
is a sample printout from the system. 

h. HTjmidity Measuring Equipment 

The Hygrodynamics Digital I hygrometer indicator 
and the Dunmore-t3rpe Lithium Chloride sensor were used to 
measure the Relative Humidity. The equipment operates on 
the basis of resistance change within an electrolytic 
solution generating a reference voltage variance which is 
proportional to the relative hvimidity change. For further 
details see Ref. 9. 
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FIGURE 6 Sample Computer Output Showing "Best Fit" Gaussian Curve 
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III. THE EXPERIMENT 

A. EXPERIMENTAL GOALS 

The original goals of the experiment were threefold. 
First, the main goal was the development of the equipment 
and procedures to measure extinction and to demonstrate 
that extinction could be measured in the ocean environment. 
Secondly, it was desirable to correlate the observed ex- 
tinction with the measured quantities of airborne particu- 
late matter and water vapor measured as humidity. Thirdly, 
it was desired to measure the extinction caused by fog and 
rain. This last goal was not achieved since, during the 
scheduled shipboard trials, the weather was persistently 
clear. 

B. EXPERIMENTAL PROCEDURE 

1. 22-23 January 1976 

On these dates the first equipment trials were 
performed in the basement hallway of Spanagel Hall, Naval 
Postgraduate School. The laser source and monitoring 
equipment were mounted on a portable equipment rack and 
positioned at distances of 15.8 m, 31.0 m, and 46.3 m from 
the receiving equipment which was kept stationary. The 
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purpose of these trials was to develop technique and 
initial equipment settings. At that time the signal 
detector was a PIN Photodiode which, for reasons to be 
explained in Section III.B.5, was replaced by the Avalanche 
Photodiode previously described. In addition to the devel- 
opment of techniques and equipment settings, these trials 
established the Laser Source Strength Monitor reference 
reading at 31 mv. 

2. 26 January 1976 

On this date the previous trials were repeated with 
similar results. In addition, the effects of "Fan-Scan" 
amplifier settings and PAR amplifier settings were measured. 
It was determined that "Fan-Scan" amplitude had a definite 
effect on the monitor pulse height and an amplifier setting 
of 1.2 was chosen as the best setting since the response 
at settings surrounding this point was relatively constant. 
The PAR settings arrived at were: (1) Low Frequency Rolloff 
100 HZ, (2) High Frequency Rolloff 300 KHZ. These settings 
provided the cleanest signal response with minimijm pulse 
shape distortion. 

3. 2 February 1976 

On this date, the first concrete attempt was made 
to determine the equipment response to transmission through 
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an atmosphere containing essentially no water droplets. 

The trial was again made in the basement hallway at trans- 
mission distances of 15.8 m, 31.0 m, 46.2 m, 61.5 m, 92.0 m 
and 122.5 m. All equipment settings were maintained as 
previously described; however, in this trial the laser 
source was left stationary and the detector and signal 
processing equipment were positioned at the prescribed 
ranges. This procedure was utilized since it facilitated 
more rapid detector alignment. The results of this trial 



culation of the extinction coefficient as a result of the 
data to be taken during the shipboard trials of the week 
of 29 March 1976. 

4. 18 March 1976 



trip, it would be advantageous to have some measurements 
at longer distances than could be achieved in the basement 
hallway. For the calibration trial on this date the laser 
source and monitoring equipment was mounted in the small 
two-wheeled trailer, and positioned at the western end of 
O' Hare Avenue (NAVY ANNEX) in close proximity to the fuel 
station so that 115 volt A.C. power could be readily 



are plotted in Figure 7. The observed value of In 




determined here was to be utilized as In 




in the cal- 



It was decided that, prior to the 29 March Acania 
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FIGURE 7 Graph of In(IR^) vs. R for 2 February 1976 
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obtained from the distribution panel located there. (See 
Figure 8) . The detector and signal processing equipment 
was located in the instrumentation bus, which was moved 
and positioned at ranges of 91.0 m, 206.0 m, and 336.0 m 
from the source. (See also Figure 8). The results of this 
trial are plotted in Figure 9. These results together with 
the results of the 2 February trials were to be utilized 
with the data collected aboard Acania during the week of 
29 March. 

5 . Week Beginning 29 March 1976 

During this week, several extinction measurements 
were recorded aboard R. V. Acania. The majority were taken 
with the ship at anchor at a range of approximately 0.8 
nautical miles off Point Pinos. Each day the ship would 
arrive on station and anchor using radar ranges. When the 
ship settled and rode to the anchor, the exact range from 
bus to ship was determined by triangulation using surveyor's 
transits. It was during this trial, on 30 March 1976, when 
an attempt was made to obtain extinction measurements at 
consistently increasing ranges, that it was determined that 
to do so would require the use of the Avalanche Detector; 
since, at ranges greater than two kilometers, the Signal- 
to-Noise Ratio obtained with the PIN Photodiode was 
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FIGURE 8 Drawing of Experimental Arrangement Naval Postgraduate School Annex 
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FIGURE 9 Graph of In(IR^) vs. R for 18 March 1976 
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insufficient. From this determination, it was decided that 
a "calibration" trial, similar to the 18 March trial, would 
have to be made using the Avalanche Photodiode. Since this 
detector produces an extremely large signal resulting in 
saturation at close ranges, it was decided to again work 
at the Navy Annex site rather than the Spanagel Hall basement. 

6. 19-20 April 1976 

On these dates the "calibration" trials were con- 
ducted utilizing the Avalanche Photodiode. The equipment 
was again configured as described previously in Paragraph 4. 
The source strength monitor at this time also produced a 
reading of 31 iitv indicating that the laser output was re- 
maining constant. The results of the trials described here 
are plotted in Figure 10. Again, the observed value of 
In I determined here was to be utilized as In |cj 

in conjunction with the data ,to be collected during the 
final shipboard trials to be conducted the following week. 
Also, during this calibration trial, anticipating the need 
to use a large aperture opening at the longer ranges, a 
trial was conducted to determine the effect on the response 
of the detector as a function of aperture size. This was 
conducted at the longest available range (336.0 m) since 
opening the aperture greater than 0.4 cm diameter at the 
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FIGURE 10 Graph of InClR^) vs. R for 19-20 April 1976 
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closer ranges resulted in saturation. The results of this 
trial are plotted in Figure 11. As indicated the response 
is linear to an aperture radius of 1.1 cm, at which point 
saturation set in. 

7. 28-29 April 1976 

On these dates, the final data collection reported 
on in this thesis was completed. Simultaneous meteorological 
data was collected at the detector site aboard ship. The 
laser source was mounted on the 18 inch telescope in the 
bus. The bus was positioned at the water's edge on Point 
Pinos. The ship was positioned north of this location at 
various ranges from 839 meters to 6748 meters. Three 
factors arose that affected the effort to duplicate the 
exact equipment conditions of the calibration trials. First, 
because of the heavy electrical load on the 3.5 kw generator 
resulting in a low line voltage, the "Fan-Scan" amplifier 
setting had to be increased to 2.5 to maintain oscillation 
of the scanning mirror. Secondly, the source intensity 
monitor indicated some laser degradation as a result of 
the low voltage condition. The maximum monitor reading 
was 22.0 mv when the trial first began, then within 30 
minutes of operation the monitor dropped to 18.5 mv and 
remained there throughout the remainder of the trials. 
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FIGURE 11 Graph of Average Intensity vs. Aperture Size 
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This was down from a constant 31.0 mv obtained during the 
calibration trials. Thirdly, because of the long ranges 
involved, it was decided that the best procedure would be 
to set the aperture size at 4.0 cm diameter and keep that 
constant setting throughout all remaining trials. 

The equipment settings were as follows: 

PAR Amplifier: Low Frequency Roll Off, 100 Hz., 

High Frequency Roll Off, 300 kHz; 
"Fan-Scan": Amplifier, 2.5, 

Cylindrical mirror micrometer, 7; 
Laser divergence micrometer, 5; 
Detector Aperture, 4.0 cm diameter. 
The general experimental procedure for these two 
days was to maneuver the ship into position north of the 
bus at the desired range. When on station the gyro stabi- 
lized detector would be aimed at the laser source. When a 
usable signal was obtained the data collection would 
commence. The Multichannel Analyzer would be allowed to 
collect a sufficient number of counts so as to produce a 

I 

nearly Gaussian distribution. This usually took from one 
to three minutes depending upon distance from the source. 
Since it was nearly impossible to hold the ship at the 
same range throughout data collection, approximately half 
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way through the collection process the ship's range from 
the source would be determined by triangulation with the 
surveyor's transits. As soon as the data collection was 
completed the ship would transit to the next desired 
position and the process would be repeated. While the ship 
was moving to the new position the data was fed into the 
HP 9810 Computer, processed, and plotted on the associated 
plotter; thus freeing the Multichannel Analyzer to collect 
data at the next position. At the closer ranges the gyro 
would "lock-on" the laser signal; however, at ranges greater 
than about 3000 meters the signal strength was not sufficient 
for "lock-on" thus the gyro operator controlled the aiming 
of the detector manually with the "Joy-Stick" control. The 
results of the data collected and the correction factors 
resulting from the equipment changes previously described 
are reported in Section IV, RESULTS AND CONCLUSIONS. 
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IV. RESULTS AND CONCLUSIONS 



A. GENERAL 

Since the major goal of the experiment, as stated 
earlier, was the determination of the extinction coefficient 
within the marine boundary layer, the actual data gathering 
revolved about the schedule of the Research Vessel Acania. 
Because of this and the discovery on 30 March of the long 
range limitations of the PIN Diode detector, the data 
breaks down logically into two distinct experimental evolu- 
tions: (1) the 30 March experiment and trials leading up 

to it, and (2) the 28-29 April experiment and the calibra- 
tion trials leading to it. 

To facilitate the determination of the extinction co- 
efficient , > considering the paucity of long-range data, 

the Natural Logarithm of the quantity (IR^) was plotted 
versus the range, R, for each data point obtained. A 
straight line was drawn through these data points utilizing 
the ordinate intercept and slope as determined from the 
Principle of Least Squares curve fitting /Ref. 1^. The 
absolute value of the slope of this curve is equal to » 
the extinction coefficient. 
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The calibration trials served to provide short-range 
data and thus determine, to a large degree, the ordinate 
intercept of the plot; whereas, the long-range data obtained 
in the desired ocean environment, from Point Pinos to Acania, 
provided the information for the slope and consequently • 

Prior to plotting the data, the following corrections 
were made where necessary: 

1. Source Intensity Correction 

Since the majority of the data was taken with a 
laser source intensity monitor reading of 31.0 mv, the 
pulse height data taken while the laser intensity monitor 
reading differed from this was corrected either up or down 
according to the ratio of the two readings. 

2 . Aperture Size Correction 

As mentioned earlier, it was necessary to utilize 
a small aperture at short ranges and a larger aperture at 
long ranges. Therefore, the data collected at the longer 
ranges was reduced by the ratio of the squares of the 
aperture radii. 

B. THE 30 MARCH EXPERIMENT 

The 30 March 1976 experiment and the calibration trials 
proceeding it served to establish the feasibility of making 
measurements at increasing ranges as well as providing the 
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initial measurement of the extinction coefficient in the 
marine environment. The data is tabulated in Appendix C. 
Figure 12 is the plotted data including the least squares 
fitted straight line. The weather was clear and the simul- 

j 

taneous meteorological observations provided an average 

i 

relative humidity reading of 0.911 over the time of the 
shipboard measurements . was 6.2 x 10"^ meters 

C. THE 28-29 APRIL EXPERIMENT 

The experiments conducted on 28 and 29 April 1976, 
together with the preceeding calibration trials, succeeded 
in providing data at even longer ranges . The data taken 
on these two dates was considered together for calculation 
purposes because of the similar weather conditions and 
equipment operation. The data is tabulated in Appendix C. 
Figure 13 is the plotted data including the least squares 
fitted straight line. The weather was again clear and the 
relative hxmiidity on 28 April was 0.776, while on 29 April 
it was 0.774 giving an average relative humidity reading 
of 0.775 over the time of the measurements . was 3.8 x 
10 meters 



46 



I 

I 

I 

I 

I 



,0 ’ 



; 

li 

-JO, 3 

• ^ 

j:f 

ii,; 

•1 

• I ' 



I 



Pt. Pinos to ACANIA 
30 March 1976 







M 

Csi 

CO 

1 

II 

Oi 

Q. 

O 

oo 



(^ynui 



cn 







ca cQ 



cj 

a. CO 



CQ 




a cQ 
en CO 



cn 




a 



47 



FIGURE 12 Graph of InClR"^) vs. R for 30 March 1976 
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FIGURE 13 Graph of ln(IR2) vs. R for 28-29 April 1976 



D. CONCLUSIONS 



It seems possible to draw three main conclusions from 
these results. 

First, it is experimentally feasible to measure laser 
beam extinction at sea over the ocean-air interface. 

Secondly, the extinction coefficient for 0.6328 micro- 
meter laser light in clear weather, through the marine 
environment is of the order of 5 x 10”^ meters 

Thirdly, humidity may be considered as a possible para- 
meter with which to predict laser beam extinction in clear 
weather . 

The latter two items, of course, require further 
investigation. 
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V. DISCUSSION OF ERRORS 



The considerations and assumptions discussed below were 
utilized in an attempt to understand how the uncertainties 
in the measured quantities affect the uncertainty in the 
determination of^^^ . 












(17) 



the more useful form 



Jn = Z c -yun, 



(18) 



was used. Then yuc was determined to be the absolute value 
of the slope of the plot of In IR'^ versus R. From (18) 



Me - J^r - 3.^. R 
R 



(19) 



From the theory of propagation of errors /Ref. 10_7 ; 






From equation (19), 



ifi 

c) C- 



- _L ^ J- 

' -Re. > jx 7?X 



and 



_ _ 7 A I . Thus, the expression for 
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the relative error in becomes 

, r i^T (^f (21) 

yt<. ^ 



where the uncertainties in C, I, and R are assumed to be 
independent. For the determination, two types of range 

error were assxamed to exist; one because of ship drift during 
the time data was being collected, and the second because of 
surveying uncertainties. Therefore, becomes 



- 

7 ^ 





( 22 ) 



where (^) is the relative error introduced by assxjming a 
20 meter ship drift, and the ( ^ ~ ^qqq ) term arises from 
assuming a 300 meter error at 6000 meters and a lesser 
error as range decreases. Furthermore, it is assumed that 
= 0.05, and a value of 5.0 x 10 ^ is utilized 

c T 

for . 

Then, utilizing equation (21), the above assumed values, 
and ranges from 1000 to 6000 meters in 1000 meter increments, 
a plot of the relative error versus range was constructed 
(Fig. 14). 

In addition to the above investigation, the least squares 
calculations used in the curve fitting process were used to 
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determine the variance in the slope of the plotted data. 
The Standard Deviations of the data for 28 April and 29 
April were found through separate calculations; the 
variance of each day's results were then computed. 

This procedure indicated a relative error in yx, of 

167o. 
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VI. SUGGESTIONS FOR FURTHER INVESTIGATION 



The following paragraphs are intended to present several 
ideas and modifications that, if implemented, would allow a 
more accurate determination of the extinction coefficient 
under the conditions experienced during this type of experi- 
ment. Furthermore, it would be possible to investigate 
the propagation of laser produced radiation at other wave- 
lengths . 

A new scanning device, recently procured, should serve 
to stabilize the "Fan-Scan" output under varying line volt- 
age conditions . 

A second demodulator unit, presently under construction, 
will provide an input to a digital voltmeter, thus facili- 
tating a more accurate and rapid determination of the laser 
source intensity. Additionally, a method of monitoring 
the signal detection equipment performance would prove 
valuable. 

The use of a Laser Rangefinder should be considered. 

This would provide more accurate and rapid source-to- 
detector range determination. 

Ship to ship, or ship to some other waterborne platform, 
operations should be considered. Since the shoreline, across 
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which the propagation path passes in the present experi- 
mental setup, is rock strewn and relatively unprotected, 
waves breaking on these rocks produce a significant amount 
of airborne water particles and vapor. This being directly 
in front of the source can only act to produce a fallacious 
result in the determination of extinction. An experiment 
conducted entirely over water would eliminate this problem. 

Further attempts should be made to measure 0.6328 micro- 
meter radiation extinction in fog, mist, and rain. Also, 
an effort should be made to measure the extinction coeffi- 
cient of other wavelength laser radiation; specifically, 

CO 2 produced, 10.6 micrometer radiation. 
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APPENDIX A 



SPHERICAL WAVE SOLUTION TO THE WAVE EQUATION 

The following development will show that a laser, which 
produces pure monochromatic radiation, can be considered to 
be a point source bf spherical waves as long as the region 
of observation of these waves is restricted to a region 
near the beam axis and at a range from the source much 
greater than the distance off the axis. 

Consider now a laser source in which the intensity 
distribution at planes normal to the propagation is 
Gaussian. Starting with the Maxwell Equations in a homo- 
geneous charge-free medium, so that V- E=0 , the 
characteristics of the Gaussian Beam are developed /Ref. 27* 

Then 7 X ^ ^ ^ J y (A-1) 



Where 



T - o— £ , 



Also 



J// 



(A-2) 



c)€ 



Taking the curl of (A-2) and substituting (A-1) for 
YkH results in 
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and 



(A-4) 



y K (VK £) - yui ^ ^ 

where Y ^ ^ x S = T ( V- ^ ) - V V* ^ 
and V * ^ - O , 



Assuming that the field quantities vary as = 

7^< 3/,z) - that is pure monochromatic radiation 

such as from a laser - the Wave Equation (A-5) becomes 

Y^£ -hk^Cr)^ = o J (A-6) 

where ^ = U>^ 6 [ I - J 

and (A-6) allows for the possible dependence of k on 
position r. At this point the development is limited to 
the case in which k^(r) is given by 

k^ ( r ) = k^ - kk 2 r^ , (A-7) 

where k 2 is some constant and k = , where A is the 

wavelength of a TEM (Transverse Electromagnetic) plane 
wave propagating in the medium. Considering the propaga- 
tion of a nearly plane wave in which the flow of energy is 
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predominantly along the Z axis and assuming a solution to 



(A-6) whose transverse dependence is on r = + y^ 

2 

then Y in (A-6) can be replaced by 

r i 4- ±J- ^ j2. 



taking E as 

-e 

equation (A-6) becomes 



— 






-Mtz 






and operating on (A-10) with (A-8) 






where 



ib’ - i± 
^ J H 



, and the assumption is 



that the transverse variation is slow enough that 



f 

Taking ^ in the form of 









substituting into (A-11) , and using (A-8) gives 



- kr^q’-2kP' - kk- 



r'^- O. 



only. 



(A-8) 



(A-9) 



(A-10) 



(A-11) 

made 



(A-12) 



(A-13) 
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If (A-13) is to hold for all r, the coefficients of the 
different powers of r must be equal to zero. Thus 



+ kq'' + kk2 = 0, (A-14) 

and . (A-15) 

^ k 

Now according to (A-7) , if the medium is homogeneous, k 2 
can be set equal to zero and (A-14) becomes 

Q^''+ kq' =0. (A-16) 



Now if the function S(z) is introduced by the relation 



Q = 




(A-17) 



It follows from (A-16) that S” = 0 and s' = a also 
S= az + b so that from (A-17) 

q(Z) = k . (A-18) 

Following Kogelnik, /Ref. 47, it will be advantageous to 
define a parameter q (z) where 




St 



Now (A-18) may be rewritten as 



(A-19) 



Equation (A-15) may now be rewritten as 



(A-20) 



(A-21) 
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where the arbitrary constant of integration is chosen as 
zero, which amounts to fixing the phase of the final solu- 
tion at zero. Substituting (A-20) and (A-22) into (A-12) 

^ Vj . (a-23) 

Taking the arbitrary constant of integration to be purely 
imaginary which leads to waves whose energy density is con- 
fined to a region near the Z axis, and re-expressing it as 
a new constant cJo 

a. -- ^ (A-24) 

— r- ^ 

where 

-- . (A-25) 

Using (A-24) and the fact that 

the first term of (A-23) becomes 



J'i ‘ ^ 'fr ' jp~h 













(A-26) 
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Substituting (A-24) into the second term of (A-23) and 
separating into real and imaginary parts 






-r 






^4^:- TJ jyi • 



(A-27) 



Defining the following parameters 



- (.) • «•// . (it,Tl , 


(A-28) 


Tf . ,1, ^ 5-77 , 


(A-29) 


(4Aj>) , 


(A-30) 


Combining (A-26) and (A-27) into (A-23) , and utilizing the 
fact that E(;^^^*) = €. results in 


i) - - r * 





Consider now the form of a spherical wave emitted by a 
point radiator placed at Z = 0 /Ref. 37. It is given by 

S' =< (--ik 



«< :k -2k ^ ) 



(A-32) 



0 









♦ 



? 9 2 

provided that X + Y <■< Z , the paraxial limitation, is 
observed. 

The last form of equation (A-32) follows where Z is 

equal to R, the radius of curvature of the spherical wave. 

? 2 2 

Since r“^ = X + y , equation (A-32) becomes 

- r 3 • (A-33) 

Now, comparing the form of equation (A-31) to that of 
(A-33) , it can be seen that (A-31) is in the form of a 
spherical wave, provided the observer is in the paraxial 
region at a distance much greater than the off-axis 
distance, and also provided that ■ >;> 1. . 
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APPENDIX B 



PROCEDURE FOR DETERMINING AVERAGE VOLTAGE, tT , 
FROM LOG VOLTAGE DISTRIBUTION 



The following development will serve to illustrate the 
validity of taking the peak value of the nearly Gaussian 
Log voltage distribution, as modified by the standard devia- 
tion of the best fit Gaussian curve, to be a measure of the 
average observed voltage . 



Let 



where 




vz 



(B-1) 




channel number in multichannel analyzer, 
voltage into the log voltage circuit. 



and 



/( = 0 ^ 2 ' 2 /, 




( XT, and known calibration 
voltages) , 



7 /' = - '0~i 



Now, define 



So that 



•r-- ^/H 



(B-2) 



(B-3) 
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From (B-3) and (B-4) , 



J-X-- 



Xt - Vl e. . 



Now 



Pc^)^v- - pL'ic)cL'?c 



where 



^TL'^)cL'^ - ^ p{,'K)cL'X 



Thus, the average voltage can be expressed as 



— ^ V~ 

^ d-d- ) 

or using (B-7) , 
from (B-6) , 

V-- -ir f> (.■>:) 

Substituting (B-5) into (B-10) for 'IT results in 



(B-4) 

(B-5) 

(B-6) 

(B-7) 

(B-8) 

(B-9) 

(B-10) 







Now assxame 



^ /f e. 



^2 0 ^- 



(B-ll) 



(B-12) 
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and replacing p(x) in (B-11) results in 



ir -- 



r - L-r- ^ 






(B-13) 



or 



V- - 



-C±Ml^Af^l2d 

'iTo!^ \ ^ A y 

(B-14) 



Expanding the exponential: 






AT%-^ 






(B-15) 



Replacing the exponent in (B-14) by (B-15) results in 






JLO~3c 









(B-16) 



Since ^ is independent of % , it can be removed 

from under the integral, and recognizing that H\€. 



is equal to C, from (B-7) and (B-12) , equation (B-16) 
becomes 



iTo €. 



-h ^ 
r X 



(B-17) 



In actual operation, p(x) was the distribution collected by 
the multichannel analyzer. In calculating the best Gaussian 
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curve fitted to the data distribution, the computer 
generated and subsequently V” . The validity 

of this result depends solely on how well the computed 
Gaussian actually fits the experimental distribution. This 
was indicated graphically by the curve drawn onto the data 
(see Figure 6) . 

Reference 8 describes the computer program in detail. 
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APPENDIX C 



DATA 



SPANAGEL HALL BASEMENT 


2 February 1976 


RANGE (meters) 


1ti(1r2) 


31 


6.89 


62 


7.15 


100 


7.15 


122 


7.31 


150 


7.20 


200 


7.15 


250 


6.95 


275 


6.80 


NAVY ANNEX 


18 March 1976 


91 


7.23 


91 


7.26 


181 


7.01 


181 


6.83 


181 


6.99 


335 


7.04 


335 


7.06 
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NAVY ANNEX 


19-20 April 1976 


RANGE 


(meters) 


In(IR^) 


91 




8.19 


91 




8.22 


91 




8.25 


91 




8.15 


206 




8.14 


206 




8.30 


336 




7.99 


336 




8.34 


336 




8.14 


PT. PINOS 


to ACANIA 


30 March 1976 


RANGE 


(Meters) 


In(IR^) 


1604 




5.42 


2384 




5.69 


3323 




5.19 


3915 




4.74 
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PT. PINOS to ACANIA 


28 April 1976 


RANGE (meters) 


In(IR^) 


1279 


7.52 


2195 


7.53 


3288 


7.42 


4420 


6.60 


4535 


6.65 


5332 


5.75 


5425 


5.69 


6748 


5.81 


PT. PINOS to ACANIA 


29 April 1976 


RANGE (meters) 


IndR^) 


839 


7.82 


2120 


7.76 


3251 


7.36 


*4156 


7.35 


5300 


6.12 


6283 


5.89 



*The following is an example of the calculations and 
corrections applied to the value of It , obtained from 
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the pulse height distribution, to arrive at the value of 
In (IR^) tabulated in this Appendix. The particular data 
point used in this example was chosen, since all correc- 

O 

tions possible were applied in the determination of In (IR"^) 



In 



J 2 ? _ , f Average voltage from pulse height distribution 

IIR J - In ^ Pre-amnlifier Gain Setting 



31 mv 

Monitor 

Reading 



Pre-amplifier Gain Setting 

0.04 .9') 

(Aperture Radius)^ j 



X 



in ^ 


, 0.538 


31 0.04 


100 


18.5 ^ 4 ^ 


In 1 


> 1557.1244 


] 


7.35 







X (4156) 
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